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The renormalization group transformation for extreme value statistics of independent, identically 
distributed variables, recently introduced to describe finite size effects, is presented here in terms 
of a partial differential equation (PDE). This yields a flow in function space and gives rise to 
the known family of Fisher-Tippett limit distributions as fixed points, together with the universal 
eigenfunctions around them. The PDE turns out to handle correctly distributions even having 
discontinuities. Remarkably, the PDE admits exact solutions in terms of eigenfunctions even farther 
from the fixed points. In particular, such are unstable manifolds emanating from and returning to 
the Gumbel fixed point, when the running eigenvalue and the perturbation strength parameter obey 
a pair of coupled ordinary differential equations. Exact renormalization trajectories corresponding 
to linear combinations of eigenfunctions can also be given, and it is shown that such are all solutions 
of the PDE. Explicit formulas for some invariant manifolds in the Frechet and Weibull cases are 
also presented. Finally, the similarity between renormalization flows for extreme value statistics 
and the central limit problem is stressed, whence follows the equivalence of the formulas for Weibull 
distributions and the moment generating function of symmetric Levy stable distributions. 
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I. INTRODUCTION 

Extreme value statistics (EVS) is a prominent statisti- 
cal problem, which attracts the attention of an increasing 
number of physicists as well as scientists from other dis- 
ciplines. In its traditional proposition EVS concerns the 
statistics of the extremal, i.e., largest or smallest value 
in a batch of random quantities, or, in general the statis- 
tics of an ordered sequence of variables near extremi- 
ties. Such extremal quantities appear and have deter- 
mining effects in a variety of areas, ranging from those in 
physics like glasses [H , interface fluctuations and random 
walks 041] j front propagation Q through engineering B 
and hydrology to seismology @, and finance (Tolllll). 
These and other applications make EVS a remarkable 
field of multidisciplinary research, where theory and data 
analysis meets practical challenges. 

The asymptotic statistics of extreme values in sets 
of independent and identically distributed (i.i.d.) vari- 
ables has been known for a long time 0, [HI, HH • Three 
main different types of limit distributions are tradition- 
ally distinguished depending on the original distribution 
(called parent distribution) of the individual variables. 
If the density function of the parent decays faster than 
any power law either at infinity or a finite upper border 
then the asymptotic distribution is the Fisher-Tippctt- 
Gumbcl (FTG) function. If the parent distribution has a 
power-law tail then the limit distribution is called Fishcr- 
Tippctt-Frechet, and if it decays as a power function 
close to an upper bound, it is the Fisher-Tippett-Weibull 
(FTW) distribution [14]]. These cases can be concate- 
nated in a single, one-parameter function form often 
called the generalized extreme value distribution [TH, [l6| . 

The question of corrections to the asymptotic distri- 
butions obviously arise in applications, since practical 
problems involve finite data sets. Knowing the magni- 
tude and shape of the finite size corrections to the limit 
distribution is thus of primary importance when analyz- 
ing EVS in empirical data. The significance of finite size 
corrections has been r ecog nized as early as in the found- 
ing paper of the field [12| ■ Then the warning was made 
that logarithmically slow convergence may seriously ham- 
per practical applications of EVS [17[ , and the first shape 
correction of the limit distribution was calculated for a 
generalized Gaussian parent [lH . The rigorous treatment 
of finite size corrections for the EVS of i.i.d. variables can 
now be found in the mathematical literature [l9l . l20j , but 
remains poorly known in the physics community and in 
the multitude of other fields of application of EVS. 

Recently a renormalization group (RG) approach was 
proposed to describe the asymptotic behavior in EVS, 
worked out in detail for i.i.d. variables in [U [22| . 
Whereas it is widely known that the concept of the RG 
has been the most efficient method for dealing with crit- 
ical phenomena [23j . its implementations have been also 
useful for other problems in probability like for sums of 
variables, i.e. the central limit problem area 0, H3]- I n 
the field of EVS, an RG method was constructed to eval- 



uate limit distributions arising in random landscape and 
interface problems @, [25[ . 

The basic idea of RG methods in statistical physics 
is to perform some coarse-graining by eliminating vari- 
ables on one scale and defining effective variables that 
describe the problem at some larger scale. Iterating this 
procedure, in each step a finite operation, leads to the de- 
scription of the large scale behavior, possibly including 
singularities, of the system. In the generic case, criti- 
cal behavior is determined by a fixed point, while near- 
critical physical systems are influenced by "eigendirec- 
tions" of the linearized RG transformation near the fixed 
point. Among other properties, the finite size corrections 
to the infinite size behavior correspond to repelling, i.e., 
relevant, eigendirections. Concerning EVS, the landmark 
study of Fisher and Tippett [l2| on the limit distribu- 
tions can be considered as a fixed point analysis of an 
appropriately defined RG operation. This simple obser- 
vation can be extended to the linear neighborhood of 
fixed points, leading to the emergence of the finite size 
correction functions [U [22[ . The difference to the results 
in the mathematical literature is that [ljl [2(| obtained 
the rate of convergence and shape corrections for fami- 
lies of initial parent distributions, while the RG study, 
a priori, produces the shape correction functions from 
an eigenvalue problem in function space, a technically 
quite simple proposition. It turns out [2l|, [22[ that the 
index of eigenvalue is in fact the exponent of the finite 
size correction, necessarily non-positive for a vanishing 
correction. Thus the finite size corrections in EVS cor- 
respond to stable perturbations about a fixed point, as 
opposed to the RG in statistical physics based on the 
elimination of variables, where size corresponds to un- 
stable directions, but in agreement with the concept as 
expounded in probability theory by 24 1. Naturally, also 
in the RG treatment connection to the parent distribu- 
tion should be established, where the RG and the existing 
mathematical treatment meet, as done in (2ll. |22|. 

Given the manifest use of the RG approach in EVS, 
in this paper we reformulate the RG study of [2l|, [22[ 
in terms of a partial differential equation (PDE) defin- 
ing continuous flows in the space of distributions. In 
that way a simple and elegant method is obtained to re- 
cover the known fixed points, lying on a fixed line, and 
the recently evaluated eigendirections. This is a most 
manageable form of the RG transformation proposed for 
EVS, possibly useful for generalizations like to higher cor- 
rections, to other-than-extremal (order) statistics, or to 
correlated variables. Interestingly, while the PDE is nat- 
urally thought as describing smooth distributions, it is 
correctly handling discontinuities. In addition, the PDE 
representation provides us with remarkable exact solu- 
tions to invariant manifolds starting out from some un- 
stable direction near the FTG fixed point. It has been 
shown before in [22j that for some initial distributions 
close to a fixed point, but differing from it in an unstable 
cigenfunction, the distribution first goes away from the 
fixed point under iteration of the RG transform, before 
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cither coming back to the same fixed point or converg- 
ing to another fixed point. Surprisingly, as it can be 
shown simply via the present PDE approach, unstable 
perturbations near the FTG fixed point conserve their 
functional form even while getting farther from the fixed 
point, with only the amplitude and the finite-size-index 
depending on the flow parameter, before the amplitude 
again becomes small during the final convergence to the 
FTG limit. In other words, the invariant manifold start- 
ing out from and returning to the FTG fixed point keeps 
being expressed in terms of a single eigenfunction, albeit 
with a changing parameter. Paths in function space other 
than the aforementioned "excursions" can also be found 
in terms of a single eigenfunction, including those hav- 
ing a finite probabilistic mass at the upper limit of the 
support. Even such distributions with discontinuities are 
described well by the PDE and are shown to converge to 
a Dirac delta, i.e., their limit distribution in EVS. These 
paths in distribution space, involving a single eigenfunc- 
tion with running eigenvalue index, can be generalized to 
families of exactly calculated RG trajectories, where one 
starts out from a continuously weighted combination of 
cigenfunctions about the FTG fixed point. Remarkably, 
one can show that such is the form of the most general 
RG trajectories in function space, irrespective whether 
they are near the fixed line or away from it. As the fi- 
nal area explored within the RG for EVS, we show that 
exact, nonperturbative, RG paths in terms of a single 
eigenfunction with running parameters can also be given 
about FTW and FTF fixed points. A seemingly para- 
doxical feature there is that a single eigenfunction not 
proper to the fixed point function can appear in those 
solutions. 

Interestingly, if one thinks of the RG flow in more 
abstract terms, without the motivation connecting it to 
EVS, it turns out that the present RG operation is one of 
the simplest that can be defined, namely raising a func- 
tion to a given power, and then rescaling it. While in 
EVS this transformation is done on the integrated dis- 
tribution function, such an operation also appears in the 
well known central limit problem (2(1 l27j of sums of ran- 
dom variables, but there such an operation is performed 
on the moment generating (or characteristic) function. 
This leads to a formal similarity between EVS and the 
central limit problem, and so we can establish a direct 
correspondence between Weibull distributions and sym- 
metric Levy-stable laws. 

The paper is organized as follows. In Sec. [IT] we moti- 
vate the concept of RG in EVS. Section Hill contains the 
main results, where IIII Al recalls basic notions including 
the definition of the RG operation. IIII Bl gives the deriva- 
tion of the PDE describing the RG flow. IIII Cl contains the 
fixed point condition, and IIII Dl presents the eigenvalue 
problem and its solutions about the fixed point. Sec- 
tion |IV] concerns the explicit, nonperturbative solutions, 
wherein IIV Al gives exactly invariant manifolds about the 
FTG fixed point, the most general form of exact RG 
paths is studied in IIVBI and IIV CI contains the nonper- 



turbative solutions around FTW and FTF fixed points. 
Finally, Sec. [V] discusses the formal similarity between 
EVS and the central limit theorem, focusing on symmet- 
ric probability distributions in the latter case. Appendix 
O discusses equivalent forms of the most general RG path 
via a direct method, without resorting to the differential 
formalism. 



II. RENORMALIZATION TRANSFORMATION 
FOR EXTREME VALUES OF I.I.D. RANDOM 
VARIABLES 

In EVS the integrated distribution fi(x) is a most use- 
ful quantity, defined for a random variable x as 

H{x) = [ p{x')dx', (II.l) 

J — oo 

where p{x) is the probability density function. It has 
the meaning of the probability of finding the variable at 
any value below x - apart from ambiguities due to Dirac 
deltas we do not treat here separately. In the simplest 
proposition of EVS, let us consider a set of N i.i.d. ran- 
dom variables Xi, i = 1, . . . , N. The probability that the 
maximum value in the set, max(xi, . . . ,.tjv), is smaller 
than a given value x is equal to the probability that all 
the variables less than x. As the variables are 

statistically independent, this yields 

Prob(max(xi, . . . , xn) < x) — fJ- (x). (H-2) 

Often /i is called the parent distribution, whence the ex- 
treme value distributions descend. 

In what follows we give a flavor of the decimation 
method leading to the RG formalism presented after- 
wards. Let us split the set of sufficiently large N random 
variables Xj into N' = N/p blocks of p random vari- 
ables each. Denoting by yj the maximum value in the 
j th block, one has 

max(si, . . . ,x N ) =max(yi, ...,y N >). (II. 3) 

The variables yj are also i.i.d. random variables, with a 
distribution p p (y) given by 

p p (y) = pP{y). (II.4) 

The above procedure, which can be further iterated, 
may thus be thought of as a RG transformation, since 
a problem involving a large number of random variables 
is transformed into a similar problem, involving a reduced 
number of variables obeying a renormalized distribution. 
Note however that this is a quite simple example of RG, 
as the variables are independent, while the RG was orig- 
inally designed, in the context of critical phenomena, to 
deal with strongly correlated variables. Yet in EVS non- 
linearities are important even in the i.i.d. case, and singu- 
larities may appear for large N, which can be dealt with 
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easily by the RG approach. In particular, one can sim- 
ply derive universal shape corrections associated to finite 
size effects that need a thorough mathematical study in 
a direct approach. Finally we make a note of that the 
parameter p is here a priori an integer number p > 2. 
In the following, we shall however consider the case of a 
continuous variable. 



III. RENORMALIZATION FLOW 

A. The RG operation 

In previous papers [U [22| , the RG transformation has 
been introduced and studied, starting out from some par- 
ent distribution and winding up near a universal limit 
function. While there the iteration parameter p was al- 
lowed to be continuous, when studying the fixed point 
and the eigenfunctions about it, that parameter was dis- 
carded soon. An alternative approach, which as we shall 
sec is well adapted to analytical studies, is to consider 
the RG flow by means of a partial differential equation 
(PDE), obtained from the discrete RG by continuation 
in p. The RG transformation introduced in [2l|, [22j had 
the form 

{R P ri(x) = p p {a(p)x + b(p)). (111.1) 

Here the scale and shift parameters a(p), b(p), to be spec- 
ified later, are included in contrast to the native formula 
pi.4jl . They aim at eliminating the degeneracy of the 
distribution developing in the large p limit and allow the 
emergence of a limit distribution. 

As iterating n times the RG transform with a scale 
factor p corresponds to an effective RG with a scale factor 
p", it is natural to introduce the variable 

s = lnp (III.2) 

in order to parameterize the RG flow. Indeed, iterations 
of the RG flow then correspond to a linear increase in 
s. The integrated distribution after p = e s iterations 
will thus be denoted as fj,(x,s). For later convenience, 
we also write the integrated distribution into a double 
exponential form, namely 

H(x, s) = e -e~ 9{x ' s) t (IH.3) 

which defines the real function g(x,s). In this way, 
/j,(x, s) is necessarily bounded between and 1 for all real 
values of g(x, s). The parent distribution can be suitably 
taken as the function at s = 

fi(x,0) = fjt(x), (III.4a) 
g(x,0)=g(x), (III.4b) 

but we will also study trajectories in function space with- 
out reference to a specific parent, when the origin of "s" 
will be set by some other convention. 



One of the major interests of the RG approach stems 
from the fact that it facilitates the study of the con- 
vergence to asymptotic distributions, including the finite 
size corrections to this limit distributions. However, from 
the practical viewpoint limit distributions are equivalent 
up to a linear rescaling x ax+b of the random variable, 
with arbitrary finite parameters a and b. In order to lift 
this ambiguity, we impose the following standardization 
conditions on the integrated distributions (two conditions 
are needed since there are two free parameters a and b in 
the above rescaling): 

p(0,s) = e- 1 , (III.5a) 
d xt i(0,s) = e-\ (m.5b) 

In terms of the function g(x 7 s) 7 these conditions trans- 
late into 

g(0,s) = (III.6a) 
d x g(0,s) = l. (111.6b) 

This convention of standardization will be kept through- 
out the paper. 

B. Partial differential equation of the flow 

1. RG transformation in the second exponent 

We now turn to the derivation of the PDE for the RG 
flow. Firstly, with the reparameterization (|III.2j) we con- 
vert the subscript and argument from p to s in the RG 
transform yielding the renormalized integrated 

distribution as 

fi(x, s) = [R s p] (x) = fi e \a(s)x + b(s)) , (III.7) 

wherein the functions a(s) and b(s) enforce the standard- 
ization conditions pil.5[) . In terms of the function g(x, s), 
the RG transform (IITL71) can be rewritten, by taking the 
double logarithm, as 

g(x,s)=g(a(s)x + b(s))-s. (111.8) 

This is a very simple operation, a linear change of variable 
in the argument of the parent (|III.4b[) and a global shift, 
wherein, however, a(s) and b(s) are determined so as to 
satisfy the nonlinear standardization conditions (|III.6|) . 

2. Driven linear PDE 

In order to make the RG transformation more explicit 
we now compute the parameters a(s) and b(s). Setting 
x = in (|HT8|) . we get from (|III.6ap that g(b(s)) = s, so 
that 

b(s)=g- 1 (s). (III.9) 
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Note that g{x) is a monotonously increasing function, 
invertible over the support of the parent distribution. To 
determine a(s), we first differentiate (|III.8|1 with respect 
to x, yielding 

d x g{x,s) = a(s)g'(a(s)x + b(s)). (111.10) 

Setting again x = and using Eq. (ini.6b|) . we get 

a(s) = l/g'(b(s)). (III.ll) 

Since we aim at describing the continuous evolution 
of g(x, s) as a function of s, we differentiate (|III.8|) with 
respect to s and obtain 

d s g(x,s) = (a(s)x + b(s))g'(a{s)x + b(s)) - 1, (111.12) 

where the dot denotes the derivative with respect to s. 
Combining equations pil.9[) and we get 



[s) = a(a), 



(111.13) 



which we substitute into Eq. (|III.12|) . Then we use 
(IIII.IOP in order to eliminate g' from Eq. (|III.12|) . and 

if we introduce as well the notation 



T(s) = 



a(s) g"(b(s)) 



a(s) 5 ' 2 (6( S ))' 
we eventually obtain the PDE 

d s g{x, s) = (1 + j(s)x)d x g(x, s) - 1. 



(III. 14) 



(111.15) 



This equation should be taken with initial condition de- 
fined by the parent distribution through (|III.4j) . and the 
parent also determines 7(s) by 14[) . We emphasize 
that the original form (|III.7jl of the RG transformation 
produces valid (monotonically increasing in x) distribu- 
tion functions from a like initial function, thus the PDE, 
derived from (|III.7[) . also must preserve monotonicity. 
This is a property maintained even if it is not obvious 
directly from the PDE. We note furthermore, that if the 
standardization condition pil.6a[) is met then, as it is 
easy to see, the solution of the PDE automatically satis- 
fies IMI.fibjl 

In (|III.15|) the coefficient function j(s) can be con- 
sidered as an "external" forcing, which, beside the ini- 
tial condition, carries information on the parent. In this 
sense, (|III.15|) is a linear, parametrically driven PDE. 



Consistently with pil.8[) . this involves a linear change of 
variable and an overall shift, both infinitesimal as 



g(x, s + ds) = gy(l + a(s)ds^jx + b(s)ds, sj ~ds, (III. 17) 

wherein the functions 5(s) and b(s) are to be specified. 
Linearizing pII.17j) with respect to ds, we get 



d s g(x 1 s) = (a(s)x + b(s))d x g(x, s) - 1. 



(111.18) 



Using the standardization conditions pil.6[) we can hence 
determine the unknown functions a(s) and b(s). Setting 
x = yields a constant b(s) = 1, whence we recover 
the previously derived PDE (|III.15[) . where in place of 
a(s) the j(s) appears, the symbol we shall use for the 
coefficient function henceforth. As a final step in this 
reasoning, by once differentiating (|III.15[) in terms of x, 
setting x = 0, and using (jIII.6b[) we find that 



j(s) = -d 2 x g(0,s) 



(111.19) 



must hold. So in this picture the PDE is the same as 
pil.isp . wherein the coefficient function 7(5) is the ini- 
tial curvature (|III.19|) of the field g(x,s). Then one can 
try to find families of solutions of the PDE, where j(s) 
is calculated from the self-consistency condition pil.191) . 
Note that the presence of 7(s) defined by (|III.19[) makes 
the PDE nonlinear, and actually nonlocal because the 
curvature at x = affects the evolution at all a?'s, but 
without external driving. 



4- Dual interpretation 



In conclusion, the RG flow is described by (|III.15|) . 
where the coefficient function 7(s) is either an external 
forcing given by the parent g{x) as described in para- 
graph IIII B 21 or defined by the initial curvature of the 
field g(x, s) itself as in lIIIB3l The two definitions pII.14|) 
and piI.19[) arc, however, identical, as one can show by 
using the standardization (|III.6|) and the form (|III.8[> for 
the RG transformation. Thus the PDE representation 
of the RG transformation has a "dual" interpretation, 
which can be thought of as a consequence of the strin- 
gent standardization requirements. The PDE (|III.15[) is 
the central result of the present paper and the source of 
the subsequent analysis. 



3. Autonomous, nonlinear PDE 

Interestingly, the same PDE pII.15j) can be alterna- 
tively interpreted as an autonomous but nonlinear flow 
equation. For this purpose we start out from an arbi- 
trary point on the flow g{x, s) and apply an infinitesimal 
rcnormalization transformation to it as 

g{x,s + ds) = [R ds g]{x,s). (111.16) 



5. Discontinuities in the integrated distribution 

An inherent property of the distribution is the extent 
of its support, i.e. the (closure of the) range in x where 
the integrated distribution is strictly monotonically in- 
creasing. Outside of the support the probability densities 
vanish. 

Firstly we mention that limited supports arise in ap- 
plications and even the known limit distributions of EVS 
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have limited support in the FTW and FTF cases [i~2| . 
We emphasize, however, that the representation (jIII.3l) 
allows us to handle the case when the support is limited 
to some range [x*(s), x*(s)] of x wherein the integrated 
distribution /i(x,s) smoothly increases from to 1. In- 
deed there the function g{x, s) increases smoothly from 
-co to 00 and no obstacle seems to arise before the ap- 
plication of the PDE (|III.15[) . 

A more problematic situation emerges when the distri- 
bution has a "point measure" , i.e. the density has a Dirac 
delta somewhere. This means a discontinuity in the in- 
tegrated distribution and thus in the g function, making 
the PDE ill-defined at first sight. Given the fact that 
later in the paper we will encounter such distributions 
here we show that such cases can be treated consistently 
within the PDE formalism as well. 

Note firstly that if the distribution n(x, 0) = /z(x) has a 
discontinuity at Xdsc then the RG transformation pil.7[) , 
for some other parameter setting s, will result in an- 
other distribution fi(x, s) exhibiting a discontinuity at 
the transformed Xdsc(s) satisfying 

Xdsc = a(s) x dsc (s) + b(s). (III. 20) 

After differentiating and using Eqs. piI.13IIII.14p we get 

idsc(s) = -1 - 7(s) x dsc {s). (III. 21) 

This equation determines the evolution of the disconti- 
nuity point in s under the original form of the RG trans- 
formation (|III.7j) . 

Now we shall consider the PDE representation of the 
RG piI.15[) in the neighborhood of the discontinuity. As- 
suming a finite discontinuity in g of step size D we have 
approximately 

g{x,s)^D6(x-x dsc (s))+E, (111.22) 

where D, E are constants and 9{x) denotes the Heaviside 
symbol. Substituting that into (|III.15jl we obtain 

-Dx dsc {s)8(x-x dsc (s)) (HI 23) 

« (1 + j(s) x) D 6(x - x dsc (s)) - 1, 

with the Dirac delta notation 6'(x) — S(x). Equating the 
coefficients of the Dirac deltas we recover Eq. pil.211) . 

The above consideration can be extended to the case 
when the discontinuity is at one of the lower or upper 
borders of the support. There the integrated distribu- 
tion jumps from zero, or to one, which corresponds by 
(MI. 31) to jumps from -co, or to 00, in the function 
g(x, s), respectively. The infinite jumps can be conceived 
as pil.221) with diverging step size D, yielding in the end 
the same condition pll.21j) . 

In conclusion we reemphasize that, if the initial dis- 
tribution contains a discontinuity, this will move just as 
prescribed by the linear change of variable in the RG 
transformation. Moreover, the singularity analysis of the 
PDE consistently yields the equation for the moving dis- 
continuity, so the PDE representation is upheld even for 



non-smooth distributions. The consistency in the evo- 
lution of the size of the jump needs further study, but 
here we concentrated on its location since this will be 
discussed on examples later in the paper. 



C. Fixed point 

Let us first look for stationary solutions g(x, s) = f(x) 
of the PDE pII.15j) independent of s. Then f(x) is in 
fact a fixed point and must satisfy the stationary version 
of (ITTIT51) 

= (l+jx)f'(x)-l, (111.24) 

where by piI.19[) we must have a constant 7 arising in 

7 = -/" (0). (HI-25) 

Hence, using the standardization condition /(0) = 0, we 
get 

f(x; 7)=/ (i+jyy 1 dy = -ln(l+ 7 x), (IH.26) 
Jo 7 

whence the fixed point distribution is 

M(*; 7 ) = e- e ~ Ux;l) = e -(i+7-)" 1/ \ (IH.27) 

This is the well known generalized extreme value distri- 
bution, obtained here as a fixed line of the RG transfor- 
mation [3, EH- The above derivation of the universal limit 
distribution family is quite brief, and was made possible 
by the PDE-reprcsentation of the RG flow. 



D. Perturbations about a fixed point 

1. RG flow for arbitrary deviations 

In order to study the behavior of the flow off the fixed 
line, we perturb a fixed point function (|III.26|) by some 
(p(x,s), not necessarily small, and then write for it the 
flow equation based on (|III.15[) . We consider two forms of 
perturbed distribution, the first choice being a modified 
argument as 

g(x,s) = f(x + 4>(x,8)). (111.28) 

For now the argument 7 has been omitted and f{x) = 
/(x; 7) used to simplify notation. It is easy to see that 

0(O,s) = O, 9*0(0, s) = (111.29) 

satisfies the standardization condition pll.6|) . Thus, us- 
ing pil.26p . we have for the 7(5) of Eq. (|III.19|) 

70) =7-^0(0,5). (111.30) 



7 



Eventually, by substituting Eq. (|III.28p into (|III.15|) and 
using (|III.26p , after obvious manipulations we obtain the 
PDE for the perturbation function <fi(x, s) 



d s (j>(x, s) = (1 + j(s)x) d x <f)(x, s) - 7 4>{x, s) 



(111.31) 



This is the form of perturbation we shall use in what 
follows. However, as an interesting side remark, we note 
that the way to introduce the perturbation is not unique, 
as soon as calculations are considered beyond linear order 
in the perturbation. Indeed, another standard possibility 
is to modify the fixed point function additively, in the 
form 



g(%, s) = f(x) + f'(x) 0(x, s), 



(111.32) 



so that to linear order in <j){x, s) it be equivalent 
to (|III.28|) . This ansatz implies again the relations 
piI.29IIII.30p . The PDE for the t hus de fined < ft(s, s) is 
then obtained by substituting Eq. (|III.32|) into (|III.15[) 

d a (j>(x, s) = (1 + j(s)x) (d x <j>(x, s) - jf'(x) (f>(x, s)) 

+ {<y{s)-j)x. (111.33) 



We emphasize that the RG flow equations (jIII.31IIII.33P 
are exact, corresponding to the definitions piI.28IIII.32l) . 
respectively. Note that the two PDE's coincide only, if 
7(s) = 7 or 7 = 0. In the latter case 

d s <j>{x, s) = (1 + j(s)x) d x <t>(x, s) + 7(a) x, (111.34) 

which we display for later usage. 

2. Linear perturbations about a fixed point 

Since the above two definitions of the perturbations 
are equivalent to linear order in (p(x, s), the correspond- 
ing PDE's piI.31IIII.33)) s hould become identical after 
linearization. From p!I.30P we see that the smallness of 
4>(x, s) implies that 'y(s) is close to 7. After linearization 
of either PDE (|III.31IIII.33|) we wind up with 

d s <p(x, s) = (1 + 735) d x </>(x, s) - 7 <f>(x, s) 

„ (111.35) 
-x %<K0, a). 

In what follows we determine the eigenfunctions emerg- 
ing from this PDE and also discuss its general solution. 

3. Eigenfunctions 



Here we are seeking the eigenfunctions of Eq. pil.351) 
defined by the property that they evolve by a purely ex- 
ponential s-dcpendcncc. As we shall sec, these are pre- 
cisely those functions (j>(x, s) in which the s and x depen- 
dence arc factorized as 



Then the standardization (|III.29|) implies 

^(0) = 0, 
<//(0) = 0. 



(III.37a) 
(III.37b) 



The factorization (|III.36|) is however not unique, as mul- 
tiplying e(s) by a constant and dividing ip(x) by the same 
constant yields the same <f>(x, s). This ambiguity can be 
eliminated by setting the scale of ip( x ) through an addi- 
tional condition; we choose 

V>"(0) = -1. (111.38) 

Then condition (|III.30j) straightforwardly gives 

e(s)=7(s)- 7 . (111.39) 

This simple relation justifies the somewhat unusual con- 
dition (|III.38|) ; furthermore, obviously e(s) must be small 
for the linearization remaining valid. Inserting (|III.36[) 
into the linearized PDE (|III.35|) we get 

e{s)ip(x) =e(s)((l + ix)tp'(x)-Tip(x)+x). (111.40) 

This can be solved only if 

e(s) 



e(s) 



7 



(111.41) 



is a constant, whence 

e(s) oc e 1 ' 3 (111.42) 

and so 

(1 + -yx)tp'(x) = (7 + Y)tp(x) - x. (111.43) 

The homogeneous part has the solution 

i/>h(x) oc (1 + 72 ;) (7+7 ' )/7 , (111.44) 

whereas a particular solution of the inhomogeneous equa- 
tion is the linear function 



ip P (x) 



1 + (7' + j)x 

+ 7) 



(111.45) 



The combination of the two so as to satisfy ^>(0) = 
necessarily leads to the solution 

^iii) = y(y + 7 ) ' (IIL46) 

where we reinserted the parameter arguments into ip. 
Note that the condition ip'(0) = is also satisfied, as 
it should. The function ^;(x;7, 7') is defined for x in 
the range 1 + jx > 0. In the case 7 = 0, the function 
ip(x;0,j') becomes 



(x,s) = e(s)ijj(x). 



(111.36) 



Mx;i) = ^(l+ix-ei' x ) 



(111.47) 
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The eigenfunction (|III.46p is the same as obtained in [2l|, 
l22l | not using PDE. There also the empirical meaning of 
the above results has been clarified: interpreting N = e s 
as the number of variables, a fixed point limit distribution 
is assumed to be reached for N — > 00 as g(\nN,x) —> 
f(x), and then 7' < is the exponent of decay in N of 
the correction to the limit distribution. Note that this 
interpretation slightly differs from the one presented in 
Sec. mi where p = e s was understood, and we kept N 
finite while performing decimation. 

4- General solution for linear perturbations 

We now turn to more general solutions of the linearized 
Eq. pil.35j) . Starting from the observation that a linear 
superposition of different eigenfunctions is also a solution, 
after some considerations we arrive at 

(j)(x,s) = (l + ~/x)x(s)+xx'(s) 

-(l + jx) X (f(x) + s), [ ■ > 

in terms of an appropriate, single argument x( x ) func- 
tion. Such a solution satisfies by construction the stan- 
dardization conditions (|III.29|> . We propose that the 
above form is the general solution of the linearized RG 
flow equation pil.35[) . 
By the specific choice 

X (s) cx e 7 ' s (111.49) 

we recover the solution (f> proportional to the factorized 
form (|III.36|) with (|IH.42|) and the eigenfunction pil.46|) . 
Choosing now x( s ) as a linear combination of such expo- 
nential functions we obtain 

X (s) = J d 7 W)e 7 \ (HI.50) 

where ^(7') is a weight function, rendering the above 
integral finite. The formula for x( s ) i s reminiscent to 
the Laplace transform of the weight function 23(7') ■ Ob- 
viously, if positive 7' indices are allowed by the weight 
function ^(7') then x( s ) w iU grow in s, so the solution 
</>(x, s) remains perturbative only for a limited range in s 
and x. Note that we do not study the question of support, 
i.e., the range in x where the linear solution corresponds 
to a valid, increasing distribution, but understand that if 
the perturbation is small then the g(x, s) of Eq. (|III.28|) 
will be increasing in some range of x about the origin. 

IV. EXACT NONPERTURBATIVE RG 
TRAJECTORIES 

An interesting phenomenon, already recognized in (22j 
is that iterating the RG transform starting from a dis- 
tribution close to the FTG distribution can lead to an 
excursion in function space quite far from the FTG dis- 
tribution, before eventually coming back to it. In other 



words, the FTG fixed point has both stable and unstable 
directions, and unstable manifolds generated by unstable 
directions loop back along stable directions towards the 
same fixed point. Concerning the FTW and FTF fixed 
points, with 7 ^ 0, a trajectory starting out from it also 
returns to a fixed point, but the latter may or may not 
be the same as the starting one. 

Remarkably, by means of the PDE formalism one can 
show that the aforementioned RG trajectory, starting 
near the FTG fixed point with an initial perturbation 
along a single unstable eigenfunction of the FTG class re- 
mains within the subspace of the eigenfunctions all along 
the trajectory, even in the non-linear regime, before it 
eventually returns to the same fixed point. Actually, this 
looping solution is one of several different families of those 
special RG paths, whose full trajectories are specified by 
a single eigenfunction with running perturbation e and 
eigenvalue parameter 7'. 

Subsequently we consider nonperturbative solutions 
obtained from a reparameterization of the general solu- 
tion of the linearized equation about the FTG function. 
It will turn out that this case in fact represents all gen- 
eral RG trajectories with arbitrary parent distribution, 
i.e., initial condition. We also clarify the relationship 
between this reparameterization and the basic RG trans- 
form, which may also be thought of as a reparameteriza- 
tion of the parent distribution. 

Finally, further trial functions for trajectories involving 
a single eigenfunction about a general fixed point with 
7^0, containing running parameters e, 7, and 7' will 
also be studied and examples given when they represent 
exact RG flows. 



A. RG flow in terms of a single eigenfunction 
about FTG 

1. The distribution and its support 

We shall now proceed to show how explicit RG trajec- 
tories emerge, where the deviation from the FTG fixed 
point can be given in terms of a single eigenfunction 
with an appropriate parameterization. That way the RG 
transformation in distribution space reduces to a two- 
dimensional flow for the perturbation strength e and the 
eigenvalue parameter 7', which we now set out to deter- 
mine. 

Let us consider the ansatz for the function g(x, s), de- 
viating from the FTG fixed point /(x; 0) = x, as 

g(x,s)=x + e(s)i>(x; 1 , (s)), (IV.l) 

where the eigenfunction ip(x; 7') is given by (|III.47j) . Here 
the function g(x,s) depends on s only through the func- 
tions e(s) and 7'(s), to be specified in the next subsec- 
tion. We remind the reader that the function g{x 1 s) de- 
termines through (|III.3[) the distribution, so if e(s) is not 
small then the distribution will depend on it nonlincarly. 
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In order to fully specify the distribution, also its sup- 
port needs to be given. The function g(x,s) as defined 
in (|III.3p was supposed to be monotonically increasing in 
x from — oo to +00 over the support, however, the above 
functi on (|TVTTj) may violate this condition. Inspection of 
(lIV.ip shows that one extremum can develop at 



1.1 7 

Zextr = — m I H 

7 



(IV.2) 



if the argument of the logarithm is positive. The x extr is a 
maximum, if e > max(— 7', 0) and it is a minimum, if e < 
min(— 7,0). In both cases there is a semi-infinite range 
in x where g{x, s) decreases, which we exclude and define 
the entire region of monotonic increase as the support. 
In particular, denoting by x* and x* the lower and upper 
border of the support, respectively, we obtain 



—00, if e > min(— 7', 0), 

£ oxtr , if e < min(-7',0), 

00, if e < max(— 7', 0), 

a; oxtr , if e > max(-7',0). 



(IV.3a) 
(IV.3b) 



Note that x ex tr depends on s through the parameters e(s) 
and j'(s). In the excluded region we simply define the 
integrated distribution n(x, s) as constant 



u(x s) = j °' if X < X *> 
W) - I 1, if x > x*, 

and, equivalently, by (|III.3[) we must have 



g(x,s) 



—00, if x < x*, 
00, if x > x* . 



(IV.4) 



(IV.5) 



By this truncation we arrive at a valid, monotonically 
increasing distribution. 

The connection to the PDE is upheld such that it 
determines the evolution of the monotonically increas- 
ing part of g(x,s). Then the evolution of the possibly 
finite extremum point :r e xtr(s) is also given, and there 
we simply understand in g{x, s) the jump from —00 to 
<?(^extr(s), s) if it is minimum, or, if it is a maximum, 
from g(x extr (s), s) to +00. In order to maintain consis- 
tency with the PDE, as it has been discussed in Sec. 
IIII B 51 such a discontinuity point has to satisfy pil.211) . 
a condition we check later. 

Finally we mention that not included in formula (|IV.3|) 
are e = 0, which is the FTG case with support being the 
real axis, and e = —7', where g(x, s) is monotonically in- 
creasing over the real axis, but does not go from —00 to 
+00 as needed by (|III.3j) to give a valid distribution. In 
other words, in the latter case the density is not normal- 
ized to unity. Here a finite limit of the support and the 
necessary jump in g{x, s) can also be introduced consis- 
tently with (|III.2ip . so that way, while there is an arbi- 
trariness in truncating the support, a valid distribution 
can be associated also to this case. 



2. Flow equations and their solution 

We know that for |e| small and 7' > the function 
(lIV.ip represents an unstable direction near the FTG 
fixed point, thus e(s) cx e S7 so long as it can be con- 
sidered small. Now we determine under what conditions 
the function form pV.ip remains valid even beyond the 
region linear in e. This ansatz is a special case of either 
pil.28j) or piL32j) . identical for 7 = 0, and with 



(x,s) = e(s)ip(x;j'(s)). 



(IV.6) 



The definition pll.30j) . together with the standardization 
PII.38[) . now specializes to 



7(s) = e(s). 



(IV.7) 



Substituting the above 4>(x, s) and 7(s) into the evolution 
equation (|III.34|) we arrive at 



e i> + e 7' dyip — (1 + ex) e d x ip 



0. 



(IV.8) 



From the expression pll.47j) of ^(2;; 7'), we obtain 

d 1 ,^ = --^+-d x ^ (IV.9a) 
d x il> = 7 V - x. (IV.9b) 
We thus end up with 



2^- 

i 



£7 



xd x ip = 0. (IV.10) 



As the functions ip and xd x ip are linearly independent, 
the expressions between brackets in the above equation 
must vanish, leading to two coupled evolution equations 
for e and 7'. After rearrangement we wind up with 



e = 2e" 

• / / 

7 = e 7 



(IV.lla) 
(IV. lib) 



Thus the RG path for the ansatz pV.l|) has been led back 
to a system of two coupled nonlinear ordinary differential 
equations. Note that on the left-hand-side we have the 
change rate of the parameters in terms of s, the logarithm 
of the scale parameter p as defined in pll.2[l . Thus the 
right-hand-side can be considered as the two-dimensional 
"beta function" for our problem, a term used in general 
for the velocity field of the evolution of coupling param- 
eters in RG methods J23|. 

While the mere fact that the ansatz pV.l[) is an explicit 
solution for RG paths implies that the discontinuity point 
evolves in s by the appropriate linear transformation of 
the coordinate, it is worth checking the consistency con- 
dition (IHI.21j) for it obtained in Sec. IIII B 51 Differenti- 
ating pV.2p by s and substituting the derivatives from 
dTvTTI) yields 



ioxtr = -^7 In ( 1 + -j } 



(IV.12) 
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Using (|IV.2p again we recover just the evolution equa- 
tion (|III.2ip of a discontinuity. Thus the ansatz (|IV.1|) 
with a discontinuity at the extremum can be consistently 
understood as representing a valid distribution with a 
discontinuity. 

Turning to the solutions of pV.lip . some of them can 
be immediately seen at a glance. Firstly, all fixed points 
lie on the e = fixed line. 

Secondly, 7' = is an invariant line, where Eq. 
pV.llap simplifies to 



2e' 



(IV.13) 



Note that Eq. pV.lOp becomes singular for 7' — > 0, but 
if 7' = is inserted in PV.8P then we again arrive at 
pV.13p . Its solution with initial condition eo = e(0) is 



e(s) 



c-o 



1 - 2e s 



(IV.14) 



Thus we obtain the remarkable feature that the FTG 
fixed point perturbed by the marginal eigendircction, i.e., 
PV.1P with 7' = 0, is a solution not only when the 
perturbation is small, but also when it is of arbitrary 
magnitude. One should keep in mind, however, that 
the support of the distribution is given by pV.3P with 
a^extr = l/ e ('S), thus in this case it does not extend to the 
full real axis. 

Thirdly, by substituting e = —7' we again arrive at an 
invariant line, where 



thus 



e(s) 



1 - e s 



(IV.15) 



(IV. 16) 



This case corresponds to the somewhat pathological sit- 
uation discussed in the end of the previous subsection, 
where the function g(x, s) is monotonic over the entire 
real axis but does not develop the divergences in both 
directions ±00 as needed to yield a valid distribution. 

In the general case, we can determine the relation 
6(7'), whose differential equation can be obtained from 
Eq. (IWTll as 



— - — + 1 

dry' 7' 



(IV.17) 



The homogeneous part of this linear differential equation 
has the general solution = A'y' 2 , while a particular 
solution of the inhomogeneous equation can be given as 



-7'. Hence in general 

e = AY 



7 



(IV.18) 



where the constant is determined by the initial condition 



£ o + 7o 
7o 2 



(IV.19) 



Now we can obtain the differential equation for 7'(s) by 
substituting pV.18|) into pV.llbp as 



j' = j' 2 (Aj'-l). 

This solves to 

1 1 , eo + 7o 



1 7o 



7o 2 



In 1 



co c -y' 



(IV.20) 



(IV.21) 



where pV.19p had been used. Hence by Eqs. 
PV.18IIV.2ip we have the sought solution in the form 
of explicitly given 6(7') and 3(7') functions. 



Trajectories in parameter space 



From Eq. (IIV. 18|) we see that the RG flow takes place 
generally on parabolas in the plane (7',e), with the ex- 
ception of the straight lines e = —7' and 7' = 0. The first 
line is obtained when A = 0, while the second one for in- 
finite \A\ with the limits A — > ±00 corresponding to the 
positive and negative parts of the ordinate, respectively. 

Typical paths in parameter space are plotted in Fig. 
Q] The direction of the flow is also marked as determined 
from pV.llbp : the 7' grows in the quarter planes I, III 
and decreases in II, IV. The flow remains in the same 
quarter-plane as it started from. The following main 
cases can be distinguished (references to the labels of var- 
ious line segments will be indicated in parentheses like (a) 
for segment "a" etc.). 

(i) Convergence in the plane - excursions in function 
space: (a) In quarter IV, with < — eo < 7 , i-e. A > 0, 
the distribution converges to FTG along the marginal 
direction 7' — > 0. The initial condition can be arbitrar- 
ily close to FTG, differing from it in an unstable eigen- 
function. (b) In quarter II, with < eo < — 7 , the 
path converges to FTG along a linearly stable direction 
7' — »• 1/A < 0. It can start arbitrarily close to FTG, 
differing from it in the marginal eigenfunction. In both 
cases (a,b) Eq. PV.3P shows that the support extends to 
the real axis. 

(ii) Convergence in the plane - non-returning solutions 
in function space: (c) In quarter IV, with — eo > 7 > 0, 
i.e. A < 0, the path converges to FTG along the marginal 
direction, (d) In quarter III the trajectory converges to 
FTG along a linearly stable direction with 7' — > 1/A < 
0. Equation PV.3P shows that in both cases (c,d), with 
the exception of — eo = 7q, i.e. A = 0, the support is 
bounded from below where the distribution exhibits a 
jump. This restriction does not affect the FTG limit for 
the statistics of the maximum, but the paths diverge in 
parameter plane when one follows them backward. From 
Eqs. pV.2IIV.3IIV.18P one sees that along the divergence 
the lower border 2* — > from below. 

(Hi) Divergence in the plane - no convergence to FTG: 
(c) In quarter II, with e > — j'o > the path di- 
verges. From PV.3P it follows that for e > i.e. 
A > 0, the support is bounded from above, with the 
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bound x* converging to zero as it can be seen from Eqs. 
pV.2|IV.3|IV.18p . There we understand a jump in the 
integrated distribution to 1. But in the case when the 
upper border of a distribution has a finite probabilistic 
weight, the limit distribution of the maximum must be 
just the discrete distribution at the limit of the upper 
border. This gives now x = with probability one. Since 
the support is unbounded from below, paths here can be 
continued backward along s — > — oo to FTG, converging 
to it along the marginal direction, (f) Starting anywhere 
in quarter I leads again to a path diverging from FTG, 
similarly to case (e) Continuing backward gives FTG, 
asymptotically along a linearly unstable direction. 



such that at s = the — e(s) attains its maximum. This 
happens at 7 = 7;'/2 where eo = — 7j/4. Then from 
(|IV.21|) we get 



1 2 1 , /V 
--^ + -lnU-l 



i 7i 7i VT 



(IV.22) 



From this equation, the asymptotes of 7'(s) can be com- 
puted as 



7i '(l-e-l^'), if a 



(IV.23) 



if s — > 00. 
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FIG. 1: Main types of RG paths when the deviation from 
FTG is given by a single eigenfunction as in (|IV.40I) . The 
characteristic cases of A = —1,0, 1 together with the limits 
A — ±00 are plotted. Circles mark their intersections with 
the fixed line e = 0. The segments a,b represent closed loops, 
c,d converge to FTG, whereas e,f diverge from it. 



In what follows we will give some explicit examples for 
the different types of RG trajectories discussed in this 
subsection. 



4- Excursions in distribution space 

We now illustrate the solution given in the previous 
subsection on the case of the returning invariant mani- 
folds discussed in (i) of the previous subsection. If we 
start an excursion near the FTG distribution in a lin- 
early unstable direction with 1 3> — eo > then the path 
will turn back to the FTG fixed point along the marginal 
direction. This corresponds to paths like (a) on Fig. [1] 
We shall prefer a parameterization of the trajectory such 
that the s — > —00 limit corresponds to 7' — > 7/ > 0, so 
by pV.18[) we have 7j = 1/A > 0, it is the initial param- 
eter, and also e — > e\ = 0. For s — > 00 the path goes to 
the origin e = 7' = 0. There remains an arbitrariness 
in choosing the origin of the parameter s: we set it now 
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FIG. 2: Amplitude e and exponent 7' of the unstable manifold 
as function of the parameter s for the case 7; = 2. 




FIG. 3: Probability densities (|IV.24[) of the unstable man- 
ifold with the parameters given in Fig. [2] during s < 0, 
i.e., 2 > 7' > 1, where the amplitude of the perturbation 
e grows. The numbers indicated are the 7' values, the ini- 
tial value 7;' = 2 corresponds to the FTG distribution, and 
7' = 1 to the farthest deviation from it measured in e. The 
parameters associated with the decreasing sequence of 7's are 
s = -00, -1.543, -1.0681, -0.7094, -0.36950, resp. 
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FIG. 4: Same as in Fig. [3] but during s > (0 < 7' < 1) 
where the amplitude of the perturbation e decreases. The 
sequence begins at 7' = 1 and goes to 7' = 0, the latter 
corresponding to the FTG distribution because there e = 0. 
The parameters associated with the decreasing sequence of 
7's are s = 0, 0.4527, 1.0903, 2.1931, 5.0986, 00, resp. 



FIG. 5: The unstable invariant manifold starting near the 
origin in parameter plane, corresponding to the path (b) in 
Fig. [T] While 7' changes from to 7f = — 2 as marked in the 
legend, the distribution diverges from then converges to the 
FTG fixed point. 



Figure [5] displays the above obtained j'(s) and e(s) 
curves for j[ = 2. We also illustrated the corresponding 
one-parameter family of probability densities as defined 
through (jlV.ip by 

p(x;i) = d x e . (IV.24) 

First we start out from the very proximity of the FTG 
distribution and move farther from it, a sequence up to 
s = is shown on Fig. [3l then for positive parameters 
s > the e turns back towards zero and eventually the 
FTG distribution is approached, as displayed on Fig. [U 
The other looping invariant manifold corresponds to 
the line segment (b) in Fig.[TJ The solution is given by the 
same formula as pv.22jl but now with 7/ replaced by a 
7f < 0. While in the parameter plane a simple symmetry 
transformation relates segment (a) with (b), the distribu- 
tions themselves are manifestly different due to the sign 
change of the perturbation. This is demonstrated on Fig. 
[5) where the functions are oppositely skewed than those 
in the excursion on Figs. 13141 

5. Non-returning paths in distribution space 

We shall illustrate paths not making excursions from 
and to FTG, as described in paragraphs (ii) and (hi) of 
subsection IIVA3I The solution (|IV.18IIV.21jl is valid, 
with semi- infinite support for e ^ 0, whence the continu- 
ous part p c (x; 7') of the probability density function can 
be reconstructed. 



First we consider a path corresponding to the segment 
(c) of Fig. [TJ where (|IV.3p gives a finite lower border 
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FIG. 6: Continuous part of the density for the trajectory 
starting in quarter IV in parameter plane and converging to 
FTG, as it corresponds to the path (c) in Fig.[TJ The integral 
of this function is less than one, except in the FTG limit 
7' = 0. Characteristic values of 7' of the distributions are 
displayed, converging to zero while s — > 00. 



x* for 7' > 0. Then the continuous part of the density 
p c {x] 7') vanishes at and has a norm smaller than one, 
as illustrated by the sequence in Fig. [6] As explained in 
Scc. lIV ATI we associate a complementary discrete proba- 
bilistic weight with x*. This weight, never exceeding 1/e 
and vanishing in the FTG limit, could be represented by 
a Dirac delta in the density, but we omit it from the fig- 
ure. Secondly, we illustrate in Fig. [7] a trajectory of the 
segment (f) of Fig. [TJ where (|IV.3j) gives a finite upper 
border x*. For s — > — 00 we have the FTG limit, and 
for increasing s the a;* decreases and reaches zero. This 
happens at a finite s, as it can be seen from PV.21[) . 
corresponding to the point at infinity along the parabola 
segment (f ) on Fig. [TJ In this limit the non-normalized 
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FIG. 7: Continuous part of the density for the trajectory 
starting in quarter I and diverging in the parameter plane, as 
it corresponds to the path (f) in Fig. [T] The integral of this 
function is less than one, except in the FTG limit, obtained 
by backtracking for s — » — oo. Characteristic values of 7' of 
the distributions are displayed, where 7' = 2 corresponds to 
FTG, as seen from Fig. |T] 



part of the distribution becomes peaked at the origin, 
and, together with the discrete weight at the upper bor- 
der of the support (not shown on Fig. [7|), form a Dirac 
delta at the origin. This single weight is the generic EVS 
limit distribution when the parent has a discrete weight 
at the upper border of the support, and formally corre- 
sponds to the 7 — > — 00 limit of the FTW fixed point 
distribution. 

The above illustrated cases are the typical non- 
returning trajectories, converging ones go to FTG, while 
diverging ones in the plane represent convergence to a 
Dirac delta. 

This ends our examples where a single eigenfunction, 
appropriately parameterized, gives exactly the RG path 
even farther away from the FTG fixed point. 



can be reabsorbed into e through a suitable redefinition 
of s, we conclude that the ansatz pV.ljl is equivalent to 
a rescaling of x together with a reparametcrization of s. 
This suggests to use the following more general ansatz, 
where the eigenfunction is now replaced by a general so- 
lution cj) of the linearized equation as 



g(x, s)=x + cf)(f3(s)x, T](s)) , 



(IV.25) 



where /3(s) and rj(s) are functions to be determined, and 
4> is defined in (|III.48[) for 7 = as 



(x, s) = x(s) + xx'{s) - x(x + s). 



(IV.26) 



Inserting the above form of g(x, s) into the PDE (jIII.151) 
yields a set of two coupled differential equations for f3(s) 
and rj(s): 



V = P, 



(IV.27a) 
(IV.27b) 



Consequently, the ansatz (|IV.25HV.26[) is an exact RG 
trajectory if the running parameters satisfy the above 
ODE's. This is remarkable if one recalls that it repre- 
sents an extension of the solution in the linear neighbor- 
hood of the FTG fixed point by only using a suitable 
parameterization. Even more remarkably, as we shall see 
it below, the present ansatz is actually the most general 
solution of the PDE for the RG flow, and the meaning of 
the parameter functions will also be revealed. 

Before going on, we should note that the pair of ODE's 
PV.27|) can be solved explicitly for the (3{rf) and s{rf) 
functions with arbitratry initial conditions. We can ob- 
tain this solution by firstly dividing the two equations to 
yield an ODE for /3(rj), which can be integrated, and 
secondly, by substituting the thus obtained /3 (77) into 
(|IV.27aj) . which then yields explicitly s(rf). Instead of 
displaying and analyzing this solution, however, we shall 
follow a different and perhaps more surprising reasoning. 



B. RG flow in terms of the general linearized 
solution around FTG 

1. Reparameterization of the general solution of the 
linearized PDE 

Given the RG trajectories in the previous subsection 
IIV Al whose functional form is characterized by a single 
eigenfunction, it is natural to ask whether the solution 
could be generalized to linear combinations of the eigen- 
functions. In other words, assuming the form of the gen- 
eral linearized solution as described in Sec. lIII D~4"l for the 
FTG case, the question is whether we can parameterize 
it appropriately such that it provides exact RG paths. 

Considering the explicit form of the FTG eigenfunction 
(MI. 471) , one sees that 7' acts as a rescaling of the variable 
x, as well as a global rescaling factor of the function ip. 
Noticing that the variations of this global rescaling factor 



2. ODE's for a general RG path 

Let us reconsider the scale and shift parameters of 
the RG transformation, i.e., a(s) and b(s) of Sec. MI Bl 
From their definition Eqs. (ML9|IILllj) we obtain a pair 
of ODE's 

b = a, (IV.28a) 
a= -a 3 g"(b), (IV.28b) 

where we remind the reader that g(x) = g(x,0) is the 
initial condition for the PDE, corresponds to the par- 
ent distribution, and meets the standardization condi- 
tions (|III.6[) . The above equations are essentially the 
same as the formerly displayed (|IIL13IIIL14|) . The initial 
condition of the scale and shift parameters is obviously 
a(0) = 1 and 6(0) = 0. 
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Conversely, it is straightforward to show that the above 
pair of ODE's, with the aforementioned initial condition, 
can be uniquely solved resulting in piI.9IIII.lip . 

We can immediately recognize the similarity between 
the pairs of equations pV.27[) and (|IV.28p . There are 
still a few steps to make before the analogy is fully speci- 
fied. First, we should set the initial conditions in (|IV.27|) 
to 77(0) = and /3(0) = 1. Then, using these values, 
we obtain the initial function g(x) in terms of x( x ) by 
substitution of (|IV.26p into pv.25p to yield 

g(x) = g(x, 0)=x + X (0) + x X '(0) - *(*). (IV.29) 

Hence g"(x) — —x"( x ) an d we see that the two sets of 
ODE's PV.27P and PV.28P are identical. Therefore, the 
correspondence 



rj(s) = b(s) = g \s) 
Pis) = a(s) = l/g'( V (s)) 



(IV.30a) 
(IV.30b) 



holds. Eventually, it remains to be clarified what the role 
of the rj and /? parameters in formula pV.25P for the RG 
flow is. While the parameters a(s) and b(s) enter as the 
respective scale and shift in the argument of (|III.8|) . at 
first glance it is not obvious that also 77 and /? in pv.251) 
define a linear transformation. We start out from pII.8P , 
substituting the function g{x) from PV.29P and s(rj) from 
the inversion of pV.30ap . This yields 

g(x, s) =g(p(s)x + r)(s)) - s 

=px + r) + x(0) + 09a! + r])x'{0) - x{P x + rfi 

-v-vx'(o)-x(o) + x(v) 

=px(l + X '(0)) + X (V) ~ x{Px + v)- (IV.31) 
Equations PV.29P and pV.30bp have the consequence 

/3(1 + X'(0)) =1 + P X '(V), (IV.32) 
which can be substituted into pV.3ip to yield 



g(x, s) =x + Pxx'(v) + x(v) ~ X(P X + v) 

=x + 4>(P x j v)- 



(IV.33) 



The last line is obtained by the definition pV.26l) . 
whereby the starting ansatz PV.25P is recovered. 

In summary, the trial function PV.25P includes the 
most general, defining formula (|III.8[) for the RG tra- 
jectories. A direct derivation of this conclusion, without 
using the differential formalism, is presented in Appendix 
IA1 One should recall that the function X (x) in (IIV.26|) . 
if small, is related to the weight in a linear combina- 
tion of eigenfunctions, as seen in Sec. MIDI In fact, the 
motivation for pV.25P was to construct a generalization 
of the solution in the linear neighborhood of the FTG 
fixed point, valid also in the non-perturbative regime, 
i.e., for not only small X (x). While the RG paths in 
terms of a single eigenfunction about FTG, as discussed 
in Sec. IIV Al could be considered as a curiosity, now we 



see that the RG transformation generally maintains the 
functional form also for linearly weighted combination of 
eigenfunctions, if appropriately parameterized, even far- 
ther from the FTG fixed point. This is to our knowledge 
a very special property for an RG transformation, unpar- 
alleled over its areas of applications in statistical physics. 



3. ODE's for the perturbation parameter 

In order to make a connection with the single eigen- 
function case, as well as to explicitly characterize the 
deviation from the FTG fixed point, a perturbation am- 
plitude e can be introduced for a general RG trajectory. 
With the notation of Sec. HVBll lct us define 



e = P 2 x"(v)- 
Hence we can rewrite pV.27bp as 

$ = ep. 



(IV.34) 



(IV.35) 



Taking the derivative of pV.34p with respect to s, we get 
using the pair of ODE's PV.27P 



e = 2e* + P*x!"{r 1 ). 



(IV.36) 



Then reexpressing r\ from pV.34p . we can write e as a 
function of e and /? only 

e= 2e 2 + Peip(e/p 2 ), (IV.37) 

where we have introduced the auxiliary function 



x"'{x"- x {y)) 



(IV.38) 



The ODE's PV.35IIV.37")) can be considered as the gen- 
eralization of pV.lip . valid for RG paths in terms of a 
single eigenfunction. The latter is recovered as a special 
case in the following way. Consider the initial conditions 
/3(0) = 1, 7'(0) = j' , and e(0) = e , then identify 7' 
with Pj' q , choose X (x) so that X "(x) = n$e loX yielding 
if(y) = 7g, thus we are led back to pV.lip . 

Note that the parameter e(s) introduced above equals 
the "effective" j(s) of (|III.14|> . and can be considered as 
a measure of the deviation from the FTG fixed point in 
the general case, not only when a deviation is represented 
by a single eigenfunction. 



4- An example 

We shall illustrate the RG paths as given in the pre- 
vious part of Sec. IIV Bl on an example for an excursion 
from the FTG fixed point. Note that returning invariant 
manifolds were, remarkably, given in closed forms in Sec. 
IIV Al The example shown here can be considered as a 
generalization of them. Consider a linear combination 
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of unstable eigenf unctions, weighted uniformly between 
7' = and a 7q > 0, corresponding to 



1 \ 3 £ ° 
7o 



7o 



3e e 



(IV.39) 



Note that the constant 7 q now is not to be confounded 
with an initial condition, rather it is the upper border of 
set of 7' indices here. 

Choosing x( x ) as given in (|IV.39p yields through 
(1IV.29[) a parent g{x) with support extending to the real 
axis, if 27g > — 3eo > 0. The scale and shift parameters 
f3 = a, 7] = b are determined through pV.30[) and the 
perturbation parameter e through pV.34[) . This is illus- 
trated on Fig. [5] whose qualitative similarity with Fig. [2] 
is apparent. While this example is about a case when 
the support is the real axis, the general formulas in this 
section also describe cases with limited support, similarly 
to the situation when the deviation from the FTG fixed 
point was given in terms of a single eigenfunction, as dis- 
cussed in Sec. IIV Al 
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FIG. 8: Excursion from FTG starting from a linear combina- 
tion of eigenfunctions. In (|IV.29[) the function (|IV.39I) with 
eo = —1/2 and 70 = 1 was taken. The running parameters 
a(s) = /3(s) and b(s) = rj(s) in ()IV.30|I . and the effective 
e(s) defined in (|IV.37|) are displayed. Note that while Fig. [2] 
characterizes single invariant manifolds, these curves param- 
eterize a continuum bundle of manifolds. The corresponding 
probability density function sequence first diverges from then 
converges to the FTG fixed point, in a manner qualitatively 
similar to the sequences displayed in Figs. 13141 



give trajectories belonging to fixed points other than the 
FTG one. We propose such a generalization in the form 
of the trial function 

g(x, s) = f(x + e(s)ip(x; 7(a), 7 '(s)) ; , (IV.40) 



where the fixed point function / is given by (|III.26[> . the 
eigenfunction ip by (|III.46|1 . and the running parameters 
e(s), 7 (s), 7 '(s), 7 (s) are to be determined from the PDE 
(IIII.lsp . We used here the notation j(s) and j(s) to dis- 
tinguish them from the "effective" 7 (s) given in (|III.19ll . 
While the 7 should be equal to 7 for a proper eigenfunc- 
tion about the fixed point with 7, we allow them to be 
different for the present trial function. 

After substitution into the PDE (|III.15|) a somewhat 
tedious but straightforward calculation shows that such 
a proposition can be valid only if 



7(s) 
7 '(a) 



7o, 

Bj(s), 



(IV.41a) 
(IV.41b) 



where 70 and B are constants. Thus the ansatz reduces 
to 

g(x, s) = f(x + e(s)fj(x; j(s), Bj(s)) ; 70) , (IV.42) 

wherein two running parameters remained. These, again 
by 15|) . can be shown to lead to an RG path if they 
satisfy the ODE's 



e = 2e 2 + e( 7o + (B — 1)7), 
7 = 7(e + 7o -7)- 



(IV.43a) 
(IV.43b) 



Consistency with the flow about the FTG fixed point 
can be checked, if we take B — > 00 while keeping 7'(s) 
111 pV.41b|l finite, and set 70 = 0, when we immediately 
recover the pair of ODE's in (jlV.lip . 

An interesting feature emerges if we keep 70 nonzero 
while B — >• 00. Then starting near the non-FTG fixed 
point /(£;7o) but using an eigenfunction belonging to 
FTG, ijj(x; 0, 7'(s)), also leads to an exact, nonperturba- 
tive RG path in function space. 

Another remarkable property of the ansatz pV.42|) is 
that it can also describe fixed points corresponding to 
7 7^ 7 o . Indeed, it is straightforward to show that 

f(x; 7 + e) = f{x + e^(x; 7 + e, -e); 7 ), (IV.44) 

which can be considered as a special structural relation 
between the fixed point functions and the eigenfunctions. 



C. Paths about the FTF and FTW fixed points in 
terms of single eigenfunctions 

1. Trial functions for the RG path 

It is a natural question to ask whether a generalization 
of the ansatz (jlV.ip in terms a single eigenfunction may 



2. Fixed points and flow lines in parameter plane 

In what follows we shall briefly review some basic prop- 
erties of the above ODE's. There are in general four fixed 
points for nonzero 70- 

(i) e* — 7* = 0, corresponding to the fixed point func- 
tion f(x;j ). 
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(ii) e* = 0,7*= 70, again with f(x; 70). 

(iii) e* = -70/2, 7* = 0, with /(x; 7o /2). 

(iv) 7* = 7o/(5+l), e* = -Bf, with f{x- 7o/(5+l)). 
Note that in (iii-iv) the fixed point functions f(x; 7) with 
7 70 arise as made possible by the special property 
(1RT441) . 

For general flow lines in the parameter plane 7, e an 
ODE is obtained from (|IV.43j) as 



£/ 2e + £ 7(B + l)-7o 
7 7 e + 7o ~ 7 



(IV.45) 



Some invariant lines can be explicitly found as follows. 

(a) e = corresponding to staying at /(x;7o). 

(b) 7 = implying a flow by e = 2e 2 + £70, and con- 
taining the fixed points (i, iii) . 

(c) e = —Bj, along which we have 7 = 7(70 — (-8+1)7), 
containing the fixed points (i,iv). 

(d) e = (B + 1)7(7/70 — 1), containing the fixed points 
(i,iii,iv). 

For graphical illustration let us restrict ourselves to 
positive 70 and B. Then all fixed points lie in the fourth 
quarter of the parameter plane with 7 > 0, e < 0, the 
(i) is unstable, (ii,iii) are saddle points, and (iv) is stable. 
The latter attracts all RG trajectories starting within the 
quarter plane with 7 > 0, e < 0. On Fig. O we display 
some characteristic flow lines. The main difference with 
respect to the RG paths about the FTG fixed point is 
that in this example we do not have excursions, rather, 
the RG transformation connects different fixed points. 




FIG. 9: RG flow in 7, e space, when the deviation from the 
FTF fixed point with 70 = 1 is given by a single eigenfunction 
as in (|IV.42[) . with B = 2. Circles mark the fixed points (i-iv) 
listed in the text, and invariant curves of the simple types (b- 
d) as well as some qualitatively drawn flow lines are indicated. 



3. Examples 

In order to highlight the above observations we give 
two families of solutions. On the one hand, the family 



g(x,s) =fix + — ? (l + 7'a;- e 7 '*) ; 7o 



(IV.46) 



corresponds to B — > 00 while B7 = 7' and so from 
PV.43P the ODE' become 



e = 2e 2 + 6(70+7'), 
7' = 7'(e + 7o)- 



(IV.47a) 
(IV.47b) 



Note that this is the generalization of pV.lip to nonzero 
70 • As the second example we display the RG path along 
the invariant line e = —Bj for which, after some rear- 
rangements, we get 



(IV.48) 



The running parameter 7(5) satisfies the ODE given un- 
der paragraph (c). These two examples demonstrate 
manifestly that cigenfunctions of some fixed point com- 
posed with another fixed point function can lead to exact 
trajectories. 



4- Paradox: can a fixed point with an eigenfunction of 
another fixed point be a solution? 

We end this part by calling the reader's attention to 
a seemingly paradoxical situation, namely, that (|IV.42|) 
contains a fixed point function with parameter 70, while 
the eigenfunction within can have, and typically has, a 
different 7/70 parameter. Then the question arises, 
how a solution near a fixed point, i.e. having a small 
e can be expressed in terms of a non-proper eigenfunc- 
tion. The resolution of that is simple: If the eigenfunc- 
tion contains parameters changing in s then it may be 
of the functional form of the eigenfunction, but the per- 
turbation in the argument corresponding to (|III.28|) does 
not separate as (|III.36|) into an s- and ^-depending func- 
tion. So the eigenfunction with changing parameters is 
not a solution of the eigenvalue problem, thus there is 
no contradiction with the known solutions. Only if the 
s-dependence of the parameter 7 in (|IV.42p is negligi- 
ble (vanishing with e) will we obtain separability as in 
(|III.36j) . The parameter 7 can change slowly only near 
fixed points in the parameter plane, among which we are 
interested now in (i) and (ii) located on the e = axis. 
But then near the fixed point (i) the eigenfunction be- 
comes -0(x; 7(s), Bj(s)) ps —x 2 /2, a valid eigenfunction 
of /(x;7o) of index 7' = 70, furthermore, near the fixed 
point (ii) the parameter 7(s) ~ 70 so again we are facing 
a proper eigenfunction (with 7' = Bjo). In summary, 
the exact solution PV.42P does not include a case where 
an improper eigenperturbation would appear. 
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V. RG FOR CENTRAL LIMIT DISTRIBUTIONS formation on the moment generating function as 



In this section we demonstrate a remarkable albeit sim- 
ple analogy between limiting behavior in EVS and central 
limit distributions. So far we have studied the RG trans- 
formation that arose from the EVS problem, as discussed 
in Secs.HTland lllll In its original form pil.lj) the transfor- 
mation consists of a linear change of variable and raising 
the distribution function to a power. It is well known, 
however, that formally a similar operation, rescaling and 
raising to power, should be performed also in the context 
of the central limit problem on the moment generating 
function (2(| [27j. In this section we briefly review this 
similarity and show that some of the so far considered RG 
fixed points actually correspond to distributions known 
as "stable", i.e., having the property that sums of ran- 
dom variables keep the same distribution up to a scale 
transformation. The main idea of the RG for distribu- 
tions falling into the Gaussian class has been expounded 
in [24| . treated as an introductory example, where the 
eigenvalue problem for analytic characteristic functions 
has also been determined. 

Below we show that the RG flow equation introduced 
in this paper for EVS leads, with a slightly different stan- 
dardization, to a PDE whose fixed point solution corre- 
sponds to central limit distributions. We restrict our- 
selves for the sake of simplicity here to real, symmetric 
moment generating functions, so the fixed point will de- 
scribe probability densities with even symmetry. We also 
evaluate the eigenfunctions and again find a close simi- 
larity to those found in EVS. 



A. RG transformation 

Consider the random variable Z defined as a rescaled 
sum of random variables, namely 



<%,s) = [it$](g) = $ e >(s) q) 



(V.4) 



Z = 



1 N 

— y% 

a-N f-f 



(V.l) 



where Zi are i.i.d. numbers each with density P{z) and 
moment generating function 



Hi) 



(V.2) 



and ajv is a suitable scaling factor ensuring a non- 
degenerate limit distribution, if such exists. We assume 
that P{z) = P{—z), which implies (z) = 0. From the 
parity of P(z), we also have that the moment generat- 
ing function $(q) is real, and satisfies $(<?) = $(—(7). 
Then the moment generating function for the variable Z 
is known to be 



$ N (q) = * N (a N q) 



(V.3) 



also real and symmetric. Allowing a continuation of N 
and introducing s = In N we suitably define a RG trans- 



An alternative interpretation of the RG would also be to 
perform a decimation operation, similarly to the proce- 
dure presented in Sec. [IT] for EVS. The function $(g, s) 
is real and even in q, and is related to the correspond- 
ing probability density by a Fourier transformation like 
in (fvT2|) . While $(0,s) = 1 follows from dV2]) . we have 
the freedom to impose a standardization condition, which 
sets the scale factor a(s), so we require 



$(l,s) = $(-l,s) 



(V.5) 



The above transformation is in clear analogy to the RG 
operation for EVS (|III.7j) . but while RG acted on the 
integrated distribution function there, here its argument 
is the moment generating function. Note the difference 
in standardization: in EVS the distribution function be- 
came 1 at the upper limit of the support and was e _1 at 
x = 0, now 1 is reached in the origin and e~ 1 at q = — 1. 
Another difference is that we do not have a monotonic- 
ity condition in q here, nonetheless, the analogy between 
pll.7|l and (|V.4|) is obvious. 

To further emphasize the similarity to the RG equation 
in EVS, we introduce the function h(q, s) by 



$(q,s) 



(V.6) 



Then h(q) = h(q, 0) corresponds to the starting moment 
generating function $(q) = $(g,0). Since h(q,s) is even 
in q, it suffices to only consider one semi-axis in q. To 
prepare for the analogy to EVS, we shall consider the 
negative semi-axis, i.e. q < 0. The RG transformation 
(|V.4I) implies 

h(q,s) = h(a(s)q) - s. (V.7) 
From $(0, s) = 1 and the standardization |V5]) follow 



h(0, s) = 00, 
h(-l,s) = 0, 



(V.8a) 
(V.8b) 



the latter setting the scale parameter by (|V.7[) as 

h(-a(s)) = 1 + s. (V.9) 

Differentiation of (|V7j) by both variables yields 

d s h(q, s) = q d(s) h'(a(s)q) - 1, (V.lOa) 
d q h(q, s) = a(s) h'(a(s)q). (VlOb) 



Hence we obtain the PDE for h as 



d s h(q,s) = q^(s)d q h(q,s) - 1, 



where 



t(«) 



d(s) 
a{s) ' 



(Vll) 



(V.12) 
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The standardization condition (|V.5[) implies that the left- 
hand-side of the (jV.lip is zero at q = — 1, whence 



7(a) 



d q h(-l, 



(V.13) 



Thus we wound up in (|V.11|) with a nonlinear PDE, 
which resembles the PDE (|III.15|) obtained for the EVS 
problem, while the analog of (fVTT3| is (|III.19|) . The dif- 
ference can be ascribed to the modified standardization. 



B. Fixed point 

Fixed point functions for the PDE above can be found 
by looking for s-independent solutions h(q, s) = f(q) 
with a constant 7(s) = 7. From (|V.11[) we obtain 



qjf'(q) 



1. 



(V.14) 



whence the solution for q < 0, satisfying the standard- 
ization /(— 1) = 0, is 



f(q;i) = -in(-q). 

7 



(V.15) 



Hence the moment generating function in the fixed point 
is by (fVT6|) 



*(g;7) = e- { ~ q) \ 



(V.16) 



This is precisely of the form of the FTW distribution 
with 7 < originally obtained by Fisher and Tippett [l2j ■ 
Now this function corresponds to the moment generating 
function of the stable distributions with even symmetry, 
usually written as 



*a(<7) 



(V.17) 



with < a < 2. By identification, we thus get a = — 1/7. 
The restriction for a is imposed to have a valid proba- 
bility density in the relation (|V.2|) . a = 2 corresponding 
to a Gaussian and a < 2 to the Levy family. Note that 
the validity of the moment generating function (namely, 
its inverse Fourier transform should be a positive func- 
tion) is in itself not ensured for any solutions of the PDE 
(|V.11I) , but the PDE preserves validity if the initial con- 
dition was appropriate. 

In conclusion, we obtained the even limit distributions 
for the sum of i.i.d. random variables in terms of their 
moment generating function, which is on the negative 
semi-axis identical to the FTW family of limit distribu- 
tions in EVS. The correspondence holds with a = — 1/7, 
where a is the Levy index and 7 the parameter of the 
FTW distribution, in the region 7 < —1/2. 



C. Eigenfunctions 

Again in analogy to the EVS problem, we can ask 
about the eigenfunctions of the RG transformation near 



a fixed point. Consider perturbing the fixed point func- 
tion (|V.15|) , whose second argument 7 we shall omit while 
assuming it fixed, as 

h(q,s) = f(q+e(s)Tp(q)) f(q)+e(s)f'(q)TP(q). (V.18) 

The symmetry of h(q, s) in q requires the function ij>(q) to 
be odd, and standardization conditions will be imposed 
as 



<K-i) = o, 
y/(-i) = - 



whereby from (|V.13|) we obtain 

7(s) =7 + e(s) 



(V.19a) 
(V.19b) 

(V.20) 



Then substitution of (|V7T8|> into the PDE (|VTT|) and 

linear expansion in e results in 

e{a)r/)(q) = e(s) (q - 7 iftq) + g 7 tf/(q)) . (V.21) 
This can be solved for nonzero e only if 



fCf) = ' 
e(s) 7 

is a constant, whence we obtain 

7 V(?) =q-ji/j{q)+qjtlj'(q). 



(V.22) 



(V.23) 



The solution meeting the conditions (|V.19j) is straight- 
forwardly obtained, for q < it is 

7, 7') = ^ (1 + . (V.24) 

Extending for all q's with odd symmetry, and using the 
Levy index a and the notation S — —0^7' we eventually 
arrive at 



Mq) = ~y (1 - \q\ s ) 



(V.25) 



Note the similarity of the eigenfunction with the one ob- 
tained for the EVS problem (|III.46|) : both are essentially 
power functions with a linear function added to meet 
standardization. In particular, the function (|V.24j) is ob- 
tained from pll.46jl in the FTW case with a suitable shift 
and rescaling. 

The present treatise of the RG for the central limit 
problem is only meant to be a demonstration of the 
close similarity with EVS. While it can have further 
consequences for convergence to the limit, possible ex- 
act solutions of the RG equation, and generalization to 
non-symmetric distributions, we relegate those to further 
studies. 
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VI. CONCLUSION AND OUTLOOK 

The RG transformation for the EVS problem can be 
considered a simple if not the simplest nonlinear RG con- 
ceivable: it consists of raising the integrated distribution 
function to a power and applying a linear change of the 
random variable. This gave rise to a linear, first order 
PDE, albeit with a parametric driving determined by the 
parent distribution, or equivalently, the initial condition. 
On the other hand, the coefficient can also be viewed as 
the curvature of the solution in the origin, which makes 
the PDE nonlinear but without an external forcing term. 
The versatility of the PDE is reflected by the fact that 
it can also account for discontinuities in the integrated 
distribution function. 

The PDE approach allowed for a quite short and 
straightforward derivation of the fixed line of limit dis- 
tributions and of the eigenfunctions. Furthermore, we 
found various exact solutions for the RG path in terms of 
fixed points and eigenfunctions, i.e., invariant manifolds, 
valid arbitrarily far from the fixed point. In general, we 
found that all solutions, extending also to the nonlinear 
region away from the fixed line, can be interpreted as 
linear combinations of eigenfunctions with possibly con- 
tinuous weights. We suggest that this peculiar structure 
of function space, unparalleled in other RG methods in 
statistical physics, is the consequence of the extreme sim- 
plicity of the RG operation. An analogous RG operation 
in terms of a PDE could also be defined for the central 
limit problem, but on the moment generating function 
rather than the integrated distribution, and it allowed 
the simple rederivation of the known, even, Levy stable 
functions as well as the calculation of eigenfunctions. 

There is a plethora of further directions worth pur- 
suing. The most obvious application of the RG theory, 
not discussed in this paper, is to calculate finite size ef- 
fects (see [HI)- While the limit distributions for the EVS 
problem of two random variables are known, having more 
intricate properties than the above studied single vari- 
able case, an RG approach there may be useful to reveal 
finite size effects. Other natural extensions to the statis- 
tics of higher order maxima and order statistics in gen- 
eral are also conceivable. Whereas there have been many 
advances thus far for the EVS of various problems of cor- 
related variables, especially in random walks, extensions 
of the RG concept again may prove to be useful. In fact, 
in a multitude of problems in the broader field of front 
propagation, like random fragmentation, extremal paths 
on trees, and randomly generated binary trees, functional 
recursions and fixed point conditions are known. Ex- 
tensions of them to appropriate RG pictures and their 
eigenvalue analysis holds the potential that finite size re- 
laxation effects can be described in terms of eigenfunc- 
tions. Another generalization of the EVS paradigm for 
i.i.d. variables exists to cases where the batch size is large 
but is itself a random variable. Here the finite size effects 
can presumably be studied again with the help of the RG 
results. Finally, the formal analogy with the central limit 



problem holds the promise of a simple approach to finite 
size effects there as well. 
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Appendix A: Relation between alternative forms of 
the RG equations 

Whereas this paper presents, and focuses on, the differ- 
ential representation of the RG transformation for EVS, 
it is insightful to perform a direct, non-differential treat- 
ment. Thus we reanalyze the relationship of the param- 
eterization (|IV.25jl with the original, defining formula of 
the RG transformation, which is essentially based on shift 
and rescaling operations. This way we again demonstrate 
that the two formulations are equivalent. We begin our 
reasoning by reconsidering the linearized RG transform 
by a direct formulation. 

So far in the paper, linear analysis of the RG has been 
based on the linearization of the PDE. An alternative ap- 
proach is to linearize dire ctly the RG transform (jIII.81) . 
which was followed in [2l|, [22| in the calculation of eigen- 
functions, a treatment we develop below. We consider a 
perturbation <fi(x,s), like in Sec. IIVB~T1 defined by 

g(x,s) = x + <f>(x,s), (A.l) 

the corresponding parent being denoted as 

g{x) = x + 4> {x) (A.2) 

Starting from (|III.8j) we get for the RG transform 

x + <ft(x, s) = a(s)x + b(s) + (f>o(a(s)x + b(s)) — s (A. 3) 

Assuming that <p = </>n n is small, one gets by linearizing 
(|HL9ll that b(s) = s - (f> (s). Similarly, (|III.13|) yields 
a(s) = b(s) = 1 — 0o(s). Inserting these expressions 
for a(s) and b(s) in the RG (|A.3j) . and linearizing the 
resulting equation yields 

(j>iin(x, s) = (f) {x + s) - x(f>' {s) - o (s). (A. 4) 

This is precisely the solution (|IV.26p of the linearized 
PDE in the FTG case, with 4>o(x) = — x( x )- Hence this 
solution can straightforwardly be obtained by linearizing 
the RG (|III.8[) . rather than by linearizing the PDE and 
then looking for the generic solution. 
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The linearized RG transform (|A.4[) can be reformulated 

as 

4>nn(x, s) = <t>o(x + s) + a(s)x + b(s), (A. 5) 

where a(s) and b(s) determine the shift by a linear func- 
tion, on </>ij n that allows it to fulfill the standardization 
conditions 

cfe n (M = 0, cW> lta (0,s) = 0, (A.6) 

yielding 

a(s) = -<f>' (s), 6(a) - -<&)(*). (A.7) 

In this way, the linearized RG clearly appears as a 
(constant-speed) translation of the variable x, followed 
by a shift by a linear function in order to restore the 
standardization conditions. 

Let us now sec how we can get the solution of the 
full RG from the linearized one in this framework. To 
this aim, we first reformulate the full RG in terms of the 
function tfi(x,s) defined in (|A.1[) . but this time without 
making any approximation. Rearranging (|A.3j) gives 

(f>(x,s) = (f>o(a(s)x+b(s)) + (a(s)-i)x+(b(s)-s), (A.8) 

as the general form of an RG trajectory. 

Now let us use 4>n n (x, s) obeying (|A.5|) and define a 
new trial function 4>(x, s) through a rescaling of x and a 
reparametrization of s, in the following way: 

4>{x,s)=(j) ]in (p{s)x,r l {s)), (A.9) 

where the auxiliary functions f3{s),rj{s) will be deter- 
mined later. Note that this trial function satisfies the 
standardization conditions (|A.6j) . We now wish to check 



whether <fi(x, s) could possibly be a solution of the full 
RG fOf . One has 

<j)(x, s) = <j) Q {Px + 77) + a{rj)Px + b(rj). (A.10) 

Comparison with (|A.8|) leads us to set the auxiliary func- 
tions as /3(s) = a(s) and r/(s) = 6(s). Then the relations 

a(s)-l = a(rj(s))/3(s), (A.lla) 
b(s)-s = b(r){s)), (A.llb) 

should also be satisfied, which is not a priori obvious. 
However, these relations define the shift of by a linear 
function, which ensures that <f> fulfills the standardiza- 
tion condition. As the trial function <f) by construction 
meets this standardization, (jA.lip should necessarily be 
satisfied. This can also be checked explicitly as follows. 
On the one hand, since 77 = b, we know that g{rf) = s, or 
from (|A.2[) equivalently 

V + Mv)=s. (A.12) 

But this equation, when we take into account (|A.7[) . turns 
out to be just (|A.llb[) . Next, taking the derivative of 
(IA.12|) with respect to s yields rj + (f/ Q {rj)fi — 1. Recalling 
(IIII.13[) we have a = b, that is, (3 = f], leading to 

P + Mv)0 = 1-, (A.13) 

which, by our using (|A.7[) . turns out to be identical to 
(Ola) . 

So we can conclude that the rather intuitively con- 
structed trial function (|A.9|) indeed incorporates the most 
general formula (|A.8jl for the RG path, just what we 
wanted to show. 
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